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We give a short introduction to the approaches currently used to describe black holes in loop 
quantum gravity. We will concentrate on the classical issues related to the modeling of black 
holes as isolated horizons, give a short discussion of their canonical quantization by using loop 
quantum gravity techniques, and a description of the combinatorial methods necessary to solve 
the counting problems involved in the computation of the entropy. 
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1. Introduction 

One of the main motivations of any quantum gravity theory is to gain an understanding of the 
quantum behavior of black holes (BH's), explain in fundamental terms the origin of their entropy 
and discuss other quantum phenomena such as Hawking radiation. Both string inspired models and 
loop quantum gravity (LQG) can be used to model quantum (or semiclassical) black holes, describe 
the microscopic degrees of freedom responsible for their entropy and explain the proportionality 
between entropy and area (the Bekenstein-Hawking law). 

In the string theory side the pioneering work of Strominger and Vafa [1] showed that, in the 
case of extremal black holes carrying certain types of charges, the statistical entropy can be obtained 
by counting microscopic BPS states corresponding to the same charges. For large values of the 
charges the horizon area of the black hole and the number of microstates are proportional and the 
constant of proportionality is, precisely, 1/4 (in geometric units) in perfect agreement with the area 
law. For other non-extremal black holes the complementarity principle of Horowitz and Polchinski 
[2] explains the proportionality between area and entropy though it does not say anything with 
regard to the proportionality constant. These days a lot of attention [3] is being payed to precision 
studies and the comparison of subleading contributions to the statistical entropy for large areas with 
the generalizations of the Bekenstein-Hawking law suggested by Wald [4, 5]. 

A somewhat parallel development has taken place in the loop quantum gravity approach to the 
study of quantum black holes. The first attempts in this direction go back to the work of Smolin 
[6] where he suggested to model quantum black holes as inner spacetime boundaries and high- 
lighted the role of Chern-Simons theory to understand the quantization of the horizon degrees of 
freedom. Other interesting suggestions by Krasnov and Rovelli followed soon after [7, 8, 9, 10]. 
A synthesis of the foregoing ideas and the development of the modeling of BH's in LQG by using 
isolated horizons appears in [11, 12]. The starting point of this approach is the Ashtekar Hamilto- 
nian formulation of spacetimes with isolated horizons as inner boundaries [13]. By quantizing the 
resulting model in a Hilbert space obtained by tensoring a Chern-Simons boundary Hilbert space 
and a standard LQG bulk Hilbert space, it is possible to address the problem of describing the sta- 
tistical BH entropy in LQG. A practical rephrasing of the combinatorial problem associated with 
the computation of the entropy was given by Domagala and Lewandowski [14]. They also pointed 
out that all the quantum area states must be taken into account to get the quantum horizon entropy 
and not only the lowest excitations, as had been stated in [12] and considered true for several years. 
Their result affected the value of the quantum of area predicted by LQG. That value was calculated 
by Meissner [15] as well as and the asymptotic behavior of the statistical entropy as a function of 
the area leading to the Bekenstein-Hawking law. 

The era of precision counting in LQG was ushered by Corichi, Diaz Polo and F.-Borja in [16, 
17]. Their main discovery was an effective quantization of the black hole entropy for small areas. 
Despite the fact that the area spectrum in LQG has a complicated structure (and, in particular, is 
not equally spaced) the statistical entropy shows a staircase structure with regular steps of constant 
area width. The first attempts to use combinatorial methods to explain this behavior of the entropy 
-in a somewhat simplified setting- appear in the works of Sahlmann [18, 19] where, in particular, 
he explored the use of generating functions. In fact, the combinatorial problems are such that it 
is possible to solve them exactly by using number theoretic methods [20], generating functions 
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[21], group theory methods, Laplace transforms and other related ideas. A detailed account of this 
approaches can be found in [22]. The very existence of the staircase behavior of the entropy and its 
eventual disappearance in the large area limit can also be studied by using the same ideas [23] . An 
important part of the present paper will be devoted to describe them. An interesting side result of 
the previous framework is the possibility of obtaining the exact partition functions (in the so called 
area canonical ensemble). By using them it is straightforward to obtain the "true" (thermodynamic) 
entropy (as opposed to the statistical or counting entropy) and its asymptotic behavior as a function 
of the area [24] . The fact that the subleading corrections of the statistical and thermodynamical 
entropies differ may be important to understand the stability of black holes. This is so because the 
overall stability of thermal systems is associated to the concavity of the entropy. 

To end this introduction we would like to mention the recent development of an explicit SU(2) 
formulation [25, 26, 27] that relies on covariant Hamiltonian methods and avoids the partial gauge 
fixing of the standard approach [11, 12]. This recent work puts in a rigorous footing some older 
proposals by Smolin [6], Krasnov [28], and Kaul and Majumdar [29, 30] where the authors sug- 
gested that the quantum black hole degrees of freedom could be described with the help of a SU(2) 
Chern-Simons theory. An important difference between these two points of view concerns the 
treatment of the quantum boundary conditions. 

The combinatorial problems that have to be solved to compute the entropy differ in some 
details but can be dealt with by the aforementioned methods. This difference percolates to the 
physics described by each model, in particular with regard to the subdominant contributions to the 
asymptotic behavior of the entropy as a function of the area. 

2. Weakly isolated horizons and canonical gravity 

2.1 Non-expanding and weakly isolated horizons 

Consider a 4-dimensional spacetime (M,g), where M is a 4-manifold and g is a Lorentzian 
metric tensor with signature (— , +, +, +). A 3-dimensional hypersurface JV C M is called a null 
surface, whenever the pullback 1 g^ 'of the spacetime metric onto JV is degenerate. This means 
that at every point of JV there exists a non-zero vector i tangent to JV and such that 

t g { ;Px h = o, 

for every vector X tangent to J/ ' . Being normal to jV the vector field i is twist free. It is generally 
true that the curves tangent to the degenerate directions of g are (unparametrized) null geodesies in 
M (see, for example, [31]). Under mild topological conditions [31] (natural from a physical point 
of view) the null 3-surface can be thought of as the world-3-surface of a 2-dimensional spacelike 
surface S each point of which travels in M at the "speed of light" in the direction orthogonal 
to S. A spacelike section S of jV can be thought of as an embedding of the abstract 2-surface S 
corresponding to a given instant of time. Suppose that the 2-surface S is a sphere and its embedding 
S is the boundary of a region of M diffeomorphic to a space-like 3-ball. Then, the null world-3- 
surface JV is a history of a 2-surface S blowing up or shrinking at the speed of light. In general, 
the 2-metric induced on S depends on the embedding S. 

'Throughout this paper, given an arbitrary covariant tensor T in an arbitrary manifold M and an arbitrary submani- 
fold A', the pullback of T to N via the embedding N ^ M will be denoted by 7W . 
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The null surfaces considered in this paper will not be generic because we will add more struc- 
ture to them. Let us start by considering the expansion of a null vector field £ tangent to jY defined 
as follows 

6 {e) := g W" b J? i8 if, (2.1) 

where g^' ah is any tensor such that 

°ab &aa! o °b'b 

On a null surface the vector field £ is defined up to local scalings £ t-t £' = ft and the expansion 
transforms as 6m i-> dig) = f6m. Hence, the value of the expansion of a null surface at a point 
x S jY is not well defined -it depends on the null normal £- unless it is zero because 

e {e) (x) = o ^e in (x) = o. 

If 6m(x) = for all x € jY , we will say that JY is a non-expanding null surface. 

For every non-expanding null surface JY, the following physically reasonable inequality 

R ah t£ h > 0, (2.2) 

involving the Ricci tensor R a b, may be assumed to hold at every x G JY and for any null vector 
£ tangent to JY at x. If this is the case, a geometric identity called the Raychaudhuri equation 
implies 2 that 

&#% ] = 0. (2.3) 

This property independent of the choice of the null vector field £ tangent to JY in the sense that 
if it holds for one choice, then it is true for all of them. According to the Einstein equations, the 
assumption (2.2) about the Ricci tensor is equivalent to 

R ah £ a £ b = %%GT ab £ a £ b > 0. 

This is the matter energy-density positivity, a physically justified assumption, in fact, as a con- 
sequence of the Raychaudhuri equation, the stronger condition T a \,£ a £ h = holds. An example 
of a non-expanding null surface for which (2.2) is true is a null plane in Minkowski spacetime, 
therefore, there is nothing necessarily exotic about this type of null surfaces. What makes a non- 
expanding null surface JY a non-expanding horizon (NEH) are, in addition to (2.2), the following 
topological and geometrical assumptions: 

• There exists a diffeomorphism S x (0, 1) — > JY , where S is a 2-sphere. 

• The intervals {x} x (0, 1) correspond, via this diffeomorphism, to null geodesies in jY . 

• Finally, each global section of S x (0, 1) corresponds to a space-like 2-surface in jY . 

In summary, a NEH can be thought of as a 2-sphere 5 which, on one hand, is blowing up 
or shrinking at the speed of light but, on the other, satisfies that its intrinsic geometry "does not 



In a non-expanding null surface, any null vector field I is twist-free and expansion-free and, hence, the Raychaud- 
huri equation implies that it is also shear- free [32]. 
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change" in the sense that it is invariant under the action of the diffeomorphisms generated by the 
null normals. 

To finally arrive at the definition of a weakly isolated horizon we will endow a non-expanding 
horizons jY with more structure, namely, we will reduce the arbitrariness in the choice of the null 
vector field £ tangent to JV by fixing a particular one (modulo constant rescalings). To this end, let 
us recall that the following identity [3 1 ] holds on every null surface 3 

V t l = k£, (2.4) 

where K : jV — > K. This means that, in particular, one could choose £ in such a way that its 
integral curves are affinely parametrized geodesies and, then, K would vanish everywhere. Instead, 
we assume a weaker condition, namely 

K = const. (2.5) 

Notice, however, that given a null surface there is still a large family of null vector fields which 
satisfy equation (2.5). 

A weakly isolated horizon (WIH) is a pair (<yK, [£]), where jV is a non-expanding horizon and 
[£] the class (modulo constant rescalings and represented by £) of null vector field tangents to JV , 
such that (2.5) holds. Notice, however, that the value of K* is not defined on the class [£] and, hence, 
it is not well define for a WIH. 

An interesting and physically relevant example of WIH is provided by the Kerr-Newman 
spacetimes (including Schwarzschild) for which the event horizon is, in fact, a non-expanding 
horizon. In each of those cases there is a Killing vector field t, such that its restriction £ to the 
horizon is null and tangent to it. This Killing vector is defined up to a constant rescaling, however, 
this ambiguity can be removed by imposing a normalization condition at the asymptotically flat 
region of the spacetime. In addition to this, there is a geometric identity, called "the Oth law of 
black hole thermodynamics" (see, for example, [33]), which implies that for such a Killing vector 
field K is constant (for this reason K is called the surface gravity). Indeed, for every null vector 
field £ tangent to a null surface we have 



v a e b 


= co a £ b 


and the following geometric identity holds: 




dK = 


= Sftco. 



Owing to the fact that £ is the restriction of a spacetime Killing vector field to a null surface, the 
right hand side to the previous equation has to be zero. Hence K is constant and then [£] defines 
on the black hole event horizon of each Kerr-Newman spacetime the structure of weakly-isolated 
horizon. 

In conclusion, the concept of a weakly-isolated horizon is a quasi-local generalization of the 
definition of a black hole event horizon. Actually, the only global assumption is the topology of the 
space-like sections of non-expanding horizons. Almost all the theorems concerning the black holes 
admit quasilocal generalizations to the type of horizons introduced in this section [34, 35, 36, 37]. 



3 It is important here that the codimension is one. 
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2.2 The symplectic 2-form in the space of spacetimes 

Given a spacetime (M,g) consider a local co-frame e 1 , with / = 1 , . . . , 4, such that g = rjjje 1 <8) e J 
in a given chart. The constant matrix r\u and its inverse will be used below to lower and raise 
indices, respectively. In our paper we will choose the non-zero components of r\u to be 

T?ll = T] 22 = -??34 = -??43 = 1- 

Given a co-frame e 1 we define the corresponding (local) frame connection, that is, the matrix (T 7 /), 
I, J = 1, ... ,4, whose entries are differential 1 -forms satisfying 

de 1 + T I J Ae J = 0, r /7 + r 7/ = 0. (2.6) 

Suppose now that e 1 is a solution of the vacuum Einstein equations, and consider a 1 -dimensional 
family of solutions e\s) parametrized by s, such that e 7 (0) = e 1 . Let us denote 



de 1 ■- — 
ds 



s=0 



e ! (s). (2.7) 



Therefore 8e', to which we will refer to as a variation of the co-frame, is a vector tangent to the 
space of solutions to the Einstein equations at the point e 1 . For every function / on this space we 
define the operator 

Sf(e') := £ fie'is)). 
as .v=o 

For example, for every value of s we have the frame connection r 7 /(,y) corresponding to the co- 
frame forms e ! (s) and, according to the previous definition, we have 



as 



s-0 



r'j(s). (2.8) 



Given two curves in the space of co-frames with tangent vectors die 1 and b\e ! at the same point e 1 , 
we define the following 3 -form in M 

yiJKL8[i(e I Ae I )A82 ] r KL := ^uK^e 1 Ae J ) Ab\T IJ - b\{e' Ae J ) A^r") . (2.9) 

This 3-form is closed as a consequence of the vacuum 4 Einstein equations, that is 

d (^e IJKL 8 {1 (e 1 A e J ) A d^I^) = (2. 10) 



because 



\ £iJK L e J A (dT KL + T K M A T ML ) = . (2.11) 



The same is true for each element of the family of 3-forms 

UuKLdp (e 1 A e J ) A b^T^ - -8 {l (e 1 A e J ) A 8 2] T U 



This is also true for a non- vanishing cosmological constant. 
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labeled by the Immirzi parameter y, namely 

d(^e I j KL d [1 (e I Ae I )A8 2] r KL - -S^e 1 Ae J ) ASjfu] = 0. (2.12) 

Therefore, if a spacetime (M,g) satisfies the vacuum Einstein equations, and we are given a family 
of 3-surfaces £ such that every pair £1 and £2 of the family define the boundary of a 4-dimensional 
submanifold of M, then the integral 

£1(81,82) := ^j[Q e ^ 5 [i( e/Ae/ ) A5 2] r * L - ^(^A^A^j (2.13) 

is independent of the choice of the surface £. Here and in the following we relax the notation in 
the usual way and denote the tangent vectors Se 1 simply by 8. 

In a more general case, relevant for black holes, we consider a region M ex of spacetime M 
bounded on one side by a causal 3-surface Jf and a family of 3-surfaces £ such that: 

• Every pair £1 and £2, defines a 4-dimensional region of the spacetime bounded by £1, £2 
and a segment of the surface JV contained between the 2-surfaces £1 n JV and £2 n Jf . 

• For every pair of tangent vectors 81,82 there exists a 2-form OC (81,82) defined on ,jV such 
that 

( ^ukl8 {1 (e 1 A e J ) A 8 2] T KL - -8 {l (e 1 A e J ) A 8 2] T U \ = da(8 1 , 82) . (2.14) 

Under these conditions the following generalization of (2.13) 

£1(81,82) := ^ I 01(81,82) (2.15) 

r f\ „ , , ,. „ _„ 1 

Y 

is independent of £. In both cases, (2.13) or (2.15), the map (81,82) 1— > £1(81,82) defines a 
bi-linear, anti- symmetric form in the tangent space at e 1 to the space of solutions of the Einstein 
equations (on M or, respectively, M ex ). By extending this procedure to arbitrary points of the 
solution space we can define the differential 2-form e 1 (->• Q., the (pre)-symplectic form on the 
space of solutions of the Einstein equations. 

This symplectic form extends naturally to the space of pairs (e 1 ,F J j) where the 1-forms F'j on 
M are independent of e 1 as configuration variables and satisfy the metricity condition Fjj = —Tjj. 
Notice that in this first order formalism the structure equation enforcing the torsion-less character 
of F'j appears as an additional field equation. 

To end this subsection we want to point out that even though the 3-form used to define Q. has 
been introduced locally it is actually well defined globally. Indeed, if we pick an atlas on M, and 
consider points that belong to two charts, the respective co-frames are related to each other by point 
dependent 50(1,3) transformations 

e" = A'/, T"j = A I K T K L (K- l ) L j + A I K d(K- l ) K J . 

The integrands in (2.13) and (2.15) are, nevertheless, invariant with respect to these and the differ- 
ential 3-form appearing in the integrand is defined globally on M. 

7 



+ ^ jJ^euK^e , Ae J )A8 2] T KL - %(e' Ae J ) a8 2] T ij 
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2.3 The symplectic form and the weakly isolated horizons 

Let us consider a foliation of M ex by Cauchy surfaces £ intersecting a (segment of a) WIH 
{^V ,[£])■ We will calculate now the surface term over £ n JV of the integral (2.15). We will do that 
in a "gauge", that we introduce by choosing a co-frame e 1 in such a way that its dual tangent frame 
satisfies the following conditions: 

GF1. The vector field dual to e 4 coincides with a null normal belonging to the class [£], i.e. e^ = £. 

GF2. The vector fields e\ and e 2 are tangent to a foliation of jY by spacelike sections. 

GF3. The foliation is preserved by the flow of £. 

We will also introduce a function v : JY — > K, constant on the leaves of the foliation, and such that 

e a (dv) a = 1, (2.16) 

and asume that the horizon jY is bounded by the sections S Vo and 5 V] , where S v denotes the leaf of 
the foliation (a section of jV) defined by a level set of v. 

The spherical topology of the horizon sections will be relevant in the following, therefore we 
need to pay attention to the global features of the fields. In particular, there is no obstruction to 
define the co-frame 1 -forms e and e globally on JV and such that the pull back e = 0. 

However, the 1 -forms e l and e 2 (that can be used to write down the metric tensor induced on each 
S v ) can not be defined globally on jY (remember that the two sphere cannot be parallelized). In this 
case it is, nonetheless, possible to use two charts in such a way that the corresponding frames are 
related by an SO (1, 3) rotation in the overlapping region. These two charts can be chosen in such 
a way that each of them contains the entire null geodesic in the segment ,jY (the integral curves of 
£). The property (2.3) in terms of our frame reads 

^g^=^(e x ®e l +e 2 ®e 2 )^ ) = 0, (2.17) 

and the metric induced on each section S v is 

g( s = ( e l (g) e l + e 2 (g> <? 2 ) (5,,) . 

Given a WIH, equation (2. 17) tells us that all the sections are naturally isometric to each other and, 
in particular, have the same area 

a = e l Ae 2 



irrespectively of the value of v. 

Finally, the connection 1-forms F ! j must satisfy (2.6). We assume that the variations 8\ and 
&i of the co-frame in (2. 13) preserve the above gauge fixing conditions. As a consequence of (2.3), 
the property (2.17) and the gauge choice introduced above, several components of the connection 
one forms T l j vanish and, hence, the left hand side of (2. 14) becomes 

/ 1 i \ MO 

{^e IJKL 8 [1 (e I Ae J )Ad 2] r KL - -8 {l (e' Ae J ) a8 2] T u \ 

1 \ MO 

8 [{ {e l Ae 2 )A8 2] T M - -8 {l (e l Ae 2 ) A8 2] F n I . (2.18) 
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Before proceeding, we will make two remarks. 

• The first is that if we restrict ourselves to WIH's such that 

1 8T 34 a = £ a 8T l2 a = 0, (2.19) 

then the boundary term of the symplectic structure in (2.15), that is, the integral along En,/, is 
identically zero. In fact, in the case of a non-expanding horizon the previous conditions can be 
satisfied by performing suitable Lorentz rotations of the co-frame. Indeed, it is enough to enlarge 
the list of the gauge fixing conditions with the addition of the conditions 

K = 0, J^e 1 = = Jzfk 2 , (2.20) 

that can be satisfied on every non-expanding horizon by a suitable choice of the null vector field 
£ (demanding it to be the tangent to an affinely parametrized null geodesic) and obtaining the co- 
frame elements e l and e 2 by Lie dragging in the direction of £ (recall equation (2.17)). 

• The second remark is that, in the classical mechanics of WIH [36, 38, 39], the first condition 
in (2.20) cannot be imposed because £ is not fixed, that is, we admit weakly isolated horizons of all 
the non-zero constant values of K. In that case it is easy to calculate that (2.14) holds with 

a(5i,&) = 28 {l (e l /\e 2 )^d 2] (K)v. (2.21) 

From this introductory example of the Mechanics of WIH, we learn two lessons. The first is 
that in this case the degrees of freedom which contribute to the symplectic form are pure gauge 
because they correspond to the freedom of rescaling a null vector field tangent to the horizon 

£ ^f£ 

by functions which preserve the condition K = const. The second lesson is that in this case the 
boundary term of the symplectic structure -the first term in (2. 15)- takes the following form 

-i- f a(5i,&) = -L5 [ia 5q(jc)v. (2.22) 

2.4 The U{\ ) Chern-Simons formulation 



Here we will discuss in more detail the gauge fixing introduced in the preceding section. Later, 
a horizon section area a will be fixed to be constant in the relevant part of the phase space, but we 
will first allow WIH of all the possible areas and coframes constrained only by the gauge fixing 
conditions GF1-GF3 listed above. 

For the assumed class of frames the right hand side of (2.18) is 

2^(5 [1 (e 1 Ae 2 )^S 2] (K)v) - -{S^Ae^A^Fnf^- (2.23) 

To better understand the meaning of the connection component Fn pulled back to jV , consider 
first the case of a frame e\ such that 

Jz?^ 1 = J^e 2 = 0. 
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Then, the corresponding connection 1-form component satisfies 

r 12fl f = o, ^(f n ) {jr) = o 

and hence {TnY can be identified with its pullback onto a slice S v , that is, the Levi-Civita 
connection 1-form of the 2-metric tensor g^ defined on any cross-section 5 5 = S v of jV by the 
orthonormal co-frame ((e 1 )^, (e 2 )^). The general form of the coframes which satisfy the condi- 
tions GF1-GF3 of Section 2.3 is 

e A = rje B , A,B = 1,2, e 3 =e\ e A = e\ (2.24) 

where r^ = r^(x) is a SO(2) matrix depending, generically, on the point x G J/ '. Due to the 
transformation law for gauge connections, we have 

(r 1 ^) = (ry^ +dh 

where 

dh = 4j(r _1 )f. 

Therefore 

d{T\f V) =d{T\) {J/) = l -R{e^e^ ] (2.25) 

where R is the Riemann scalar curvature of the 2-geometry of S. 

It turns out that the second term in (2.23) can be written in the form corresponding to a U(l) 
Chern-Simons theory after a suitable choice of variables. Specifically, the suitable variable is a 
SO(2) connection AoxvJ/ with constant curvature equal to the average curvature of (r 1 ^) .In 
other words, such that 

dA = Ur) (e l Ae 2 )^ (2.26) 

where 

(R) = - f R e l Ae 2 . (2.27) 



aJs v 
That is, 

A = {TnV ] +w 

where w is a globally 6 defined 1-form on jV satisfying 

dw = -{R-{R)){e x Ae 2 )^\ J^w = 0„ £ a w a = 0. 

To find w we must notice first that the global existence of w^ on each S follows from the vanishing 
of the integral of the right hand side of the the above equality, namely 



l 2 J s (R-(R))e l Ae 2 = 0. 



5 Remember that any two cross-sections S v and S v ,< are naturally isometric to each other. 

6 As opposed to (r^) which is defined only locally, i.e. in the chart where the local frame (e l ,e 2 ) is defined. 
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Next, w is defined on jY by requiring it to be orthogonal to, and Lie dragged by, I. The connection 
A is defined on the same bundle as r ! 2, that is, on the orthonormal frame bundle. In order to agree 
with [40,41] we define 



V 



1 



(2.28) 



which is a £/(l) connection on the square root of the orthonormal frame bundle over S. 
Finally, in terms of the connection V it is easy to see that 

dV 



5 {l (e l Ae 2 )^ Ad 2] (F n )^ 
The average curvature (R) can be obtained by using the topological invariant 



{ W) ] 



Re 1 Ae 2 = 87T 



as 



<*> 



87T 

a 



Therefore, for variations 5i and 52 such that 

8\a = ^a = 0, 

we have 

5 {l (e l Ae 2 ) { ' r) AS 2] (r n )(' /r '> = d (£-SiV A^y) . (2.29) 

The final formulation of the theory is as follows. We consider the space of spacetimes each of 
which is the exterior of a weakly isolated horizon. These spacetimes satisfy the Einstein vacuum 
equations. The areas of all the horizons are equal to an arbitrarily fixed number a. The space- 
time geometries are defined in a region M ex of a 4-dimensional manifold M contained between 
3-surfaces JV, £q and Ej (see Figure 1). 




/ 



Figure 1: Schematic picture of the region M c: 



We will build now a coordinate system (x ,x ,v) on j¥ by pulling back coordinates on 5 x 



(vo, vi). In order to do this let us fix a diffeomorphism 

,JV -► 5 2 x(v ,V!), 



(2.30) 



11 
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where S 2 is a 2-sphere, on which we have selected a chart with coordinates (x l ,x 2 ) and v S (vo,Vi). 
The interpretation of v is clear now, M ex is foliated by 3-manifolds £ whose intersections with Jf , 
given by the condition v = constant, are diffeomorphic to a 2-sphere. We restrict ourselves to metric 
tensors g on M such that (jV , [d v ]) is a weakly isolated horizon in which the lines x 1 ' =constant 
are null geodesies affinely parametrized by v (the parametrization of [vo,vi]) and the surfaces £ 
are spacelike. Each of the allowed metric tensors g is represented by coframes e 1 defined on open 
domains covering M such that 

g = e l ®e l + e 2 ®e 2 -e 3 ®e A - e 4 (g>e 3 . 

Furthermore, each coframe satisfies at JV the conditions introduced in Section 2.3, that is, 

e\ = d v , 

e\ and e2 are tangent to the 2-surfaces 5 V defined by the condition v =const. 

According to (2.15), (2.18), (2.23) and (2.29), the symplectic form in the space of geometries 
is now 

^(<5l,<>2) =n H or(<5l,C>2) + ^Bul(<5l,52) 

where 

«W(ft,«) := 2i2WG- r l S ' VAS2V - 

We can simplify the expression of the symplectic form in the bulk (that is in E) by introducing a 
condition restricting the Lorentz transformations of the frame. This is done by assuming that e\ , e2 
and 

ejr 1 = —=(e 4 -e 3 ) 

are tangent to £, that is fixing the frame normal to £ to be 

em = -j={e4 + e 3 ). 
The bulk part of £1 takes then the following form 

^l(\ £ ^ 5 l^ lAeJ ) A5 2] rKL ~ ^(^A^A^ffj = 2 J^PfdqAi (2.31) 

where we have used real Ashtekar variables which emerge from this symplectic form quite naturally 
because 



\6nGy 



Pi = T^^. eabCe iJKin\e J )Z\e K )c\ (2.32) 



K = (r KL )P {\e l j KL + y8})n\ (2.33) 



12 



Quantum isolated horizons and black hole entropy J. Fernando Barbero G. 

where Pf stands for P", P%, and Pfr 1 and similarly for A' a . 

In this approach, the fields V on S and P, A on Z are subject to the horizon constraint [13] 



dV (S) + (4ft) G y eabc pa I £, dx b Adx cj(S) = Q (234) 



The constraint (2.34) has a clear geometric meaning. It generates the gauge transformations defined 
by the rotations at each tangent space to S 

e\ = cos2Aei + sin2Ae 2 , e' 2 =— sin2Aei + cos2Ae 2 , e'jr 1 = eir 1 (2.35) 

where A : S — > M is an arbitrary differentiable function, namely 

V = V + dA, A"n = A l n + 2JA, 

A' 1 = COS2AA 1 - sin2AA 2 , A 2 = sin2AA 1 + cos2AA 2 , 
P[ = cos2APi + sin2AP 2 , P 2 =-sin2APi + cos2AP 2 , P 7 V = Pir 1 . (2.36) 

This constraint will be used below in the integrated form, namely, for every non-self intersecting 
loop ds in S which has the interior s C S, the integrated version of the horizon constraint reads 

nl . \ ln [ V +P s , r = (2.37) 

where 

P s>r := fplr L -e abc dx b /\dx c (2.38) 

is the flux of the vector density Pj'r 1 along the 2-surface s. 

2.5 The SU (2) Chern-Simons formulation 

In addition to the U(l) framework described before there is also a recent SU(2) approach for 
vacuum gravity which admits as symmetries all the rotations on the horizon [25, 26, 27]. In this 
approach one assumes the spherical symmetry of the spacetime geometry on the isolated horizon 
and consider only horizons with a fixed area. An isolated horizon (IH) is a WIH satisfying the 
additional condition that 

[J^,D]=0, 

where D is the intrinsic derivative operator defined in jV by X"D a Y b := X"V a Y h , where X, Y are 
vector fields tangent to JV and V the Levi-Civita connection associated with g. Notice that D is 
well defined because it preserves the tangent bundle of ,jV as a consequence of ^tg^ > = 0. 

In this case the frame (gauge) fixing conditions are somewhat relaxed because neither e^ is 
assumed to be I nor e\ and e 2 are tangent to 5. However, the time gauge is still fixed, that is, n T ej 
defined in the previous subsection is still assumed to be orthogonal to Z. Upon those assumptions 
the horizon part of the symplectic form becomes 



«Hor(5i , &i) = ^-^ 2^ / 5iA' A b\A' , (2.39) 
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which corresponds to that of a SU(2) Chern-Simons theory. Finally the form of the horizon con- 
straint becomes now 



1 ( *\\ 

(-P« £abc dx a Adx b + 87r2(i a _ y2) / =°> ( 2 - 40 ) 

where Pj and A' are the same as above, and 

F l = dA'" + -e} t A-' A A*. 

The constraint (2.40) generates SU(2) gauge transformations on 5 and coincides, modulo some 
factors involving the Immirzi parameter, with the constraint given in [11] before performing the 
customary partial gauge fixing leading to the U(l) model. Notice, however [45], that the remaining 
two conditions obtained by projecting (2.40) with respect to two independent internal directions 
orthogonal to r 1 , together with the gauge fixing conditions GF1-GF3 given in Section 2.3, should be 
solved. This must be done classically (i.e. before quantizing) though this was not the path followed 
in the original papers on the subject where, instead, the additional conditions were implemented 
weakly in the quantized version. This, together with the fact that the £/(l) model does not use 
spherical symmetry, are some of the reasons that explain the discrepancy of both approaches. 

3. Quantum isolated horizons 

The Hamiltonian analysis of general relativity in the presence of isolated horizon inner bound- 
aries leads to the introduction of a Hilbert space built as the tensor product of a bulk Hilbert space 
and a horizon Hilbert space, Jff^m = ^Hor ® =^Bui- Here the bulk degrees of freedom are associated 
with the Cauchy surface £ and the horizon degrees of freedom with the intersection of £ with the 
isolated horizon S = £ n JV . The bulk Hilbert space admits a basis spanned by spin networks that 
are allowed to intersect the horizon S. These intersections are known as punctures and carry quan- 
tum numbers jj, labeling SU{2) irreducible representations associated with the intersecting edges. 
They also carry quantum numbers nij defined by projecting the spin vectors with respect to any 
space-like vector field on the horizon (that can be defined with the help of extra privileged struc- 
tures when available). In any case the horizon Hilbert space Jtf£ OT corresponds to a Chern-Simons 
(CS) theory with level k € N. As mentioned above, this may have either £/(l) or St/ (2) as gauge 
group depending on the details of the treatment. The difference between them is the different role 
played by symmetry requirements and the treatment of the quantum boundary conditions. These 
are the quantum counterparts of equations (2.37) and (2.40), respectively. They generate the gauge 
transformations on the horizon. The classical isolated horizon boundary conditions imply that the 
connection is reducible on the horizon. This means that there exist internal vectors on the sphere 
such that their covariant derivatives are zero. As explained in section 2.4, a particular choice r 1 
for this vector can be interpreted as a partial gauge fixing condition for the SU(2) connection that 
implements a symmetry reduction from St/ (2) to t/(l). This is the path followed in the original 
treatment [11], where the only component of the boundary conditions that is promoted to a quan- 
tum operator and implemented a la Dirac is the r'-projection of equation (2.40) (see [42] for a 
discussion on this issue). On the other hand, in the St/ (2) treatment the full set of conditions are 
promoted to quantum operators and enforced on the physical states without any gauge fixing. 
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In the t/(l) case [11, 12] the horizon Hilbert space J^ OI depends on the level k G N of the 
quantum CS-theory. This fixes the prequantized value of the area of the isolated horizon to the 
value au = 4ny£pk. The Hilbert space is spanned by t/(l)-CS basis states |(ci,. . . ,cjv))^ or defined 
on the punctured sphere. They are labeled by ordered sequences of non-zero congruence classes of 
integers modulo k. Each of the labels cj is an integer number in the set {1,2, ... ,k— 1} and labels 
the quantized deficit angle 4ncj/k of the 7-th puncture. The spherical topology of the horizon S 
imposes an additional restriction on the curvature that translates into the following condition 

^c/ = (mod/:), 

for the labels in a given sequence (ci, . . . ,Cff). The CS-labels (cj) are related, via the quantized 
isolated horizon boundary condition (2.37), to the quantum geometric labels mj corresponding to 
the jj representation of the edge piercing the horizon at the corresponding puncture. The condition 
that they must satisfy is 

cj = —2mj (mod k) , 

where mj G {— ji,— jj-\- l,...,jj}. This restriction on the form of the basis states |(c/))h oi .® 
\(jj,mi),-- -)bu1 °f ^Ki n is the quantum counterpart of the isolated horizon boundary condition. 

In the St/(2) proposal [25, 26] the horizon Hilbert space is that of a St/ (2) Chern Simons 
theory with a level k corresponding to the same prequantized value of the area. The quantum 
states are labeled now by representations 57 of the quantum group SU(2) q . The quantum matching 
conditions, (2.40) in this case, leads to an identification of these labels with the jj, that is, the labels 
associated with the edges of the spin network piercing the horizon at the punctures. 

4. Black hole entropy 

The quantum states of the horizon belonging to J^ or and compatible with a given value of the 
area are responsible for the black hole entropy. Given a bulk state vector labeled by a spin network 
piercing the isolated horizon, it is possible to assign an area to the horizon as the eigenvalue of the 
area operator given by 

a LQG C//,m 7 ) = 87T 7 4£ y/ji(ji+ 1) , (4-D 

/ 

where the jj are the labels of the edges at the punctures. At this point it is important to mention 
that two types of areas related to the horizon. The first is the prequantized area a\ that must be 
introduced in the quantization of the CS theory. The second is the area eigenvalue assigned to the 
horizon by the spin network and given by (4.1). In principle it would be desirable to take both 
areas as equal. However, the fact that a\ is not an element of the area spectrum precludes us from 
doing this. 7 In practice this difficulty is sidestepped by introducing an area interval \a\ — S,cik] 
with a large enough 8 ensuring the presence of a prequantized value of the area in it. Though this 
may seem an ad hoc way to solve this problem, the introduction of certain intervals is, in fact, 
customarily used in statistical mechanics to define some relevant ensembles (the microcanonical 



7 We want to point out, nonetheless, the existence of other choices for the horizon area operator within the LQG 
approach [43], with evenly spaced area eigenvalues, in which this problem disappears. 
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in particular). The ultimate reason why this is acceptable relies on the fact that, in an appropriate 
thermodynamic limit, the value of the entropy is independent of the width of the interval. Actually, 
it is possible to consider intervals of the form [ao^/t] (with «o being the lowest area eigenvalue) 
to simplify the computations. The consideration of the thermodynamic limit is also important for 
conceptual reasons related to the smoothness properties of the entropy. To be usable in standard 
thermodynamics the entropy must satisfy some regularity requirements, for example it must be 
a differentiable function of the energy (otherwise it is impossible to define the temperature). As 
shown in the classic paper by Griffiths [44] this is guaranteed in the thermodynamic limit. 

An important comment to be made at this point concerns the use of area ensembles. The 
standard framework of statistical mechanics is based on the use of energy ensembles. Both the 
standard microcanonical and canonical ensembles are defined in terms of the energy. The statistical 
entropy of a gas is obtained, for example, by counting the number of energy eigenstates below a 
given fixed energy value E. The partition function, in its stead, is obtained by adding e~P E " for 
all the possible energy eigenvalues (taking into account their degeneracies). In the case of black 
holes in LQG the role of the energy is played by the area [7]. This is somehow a necessity in the 
formalism as this is the only geometric/physical quantity that can be assigned in a natural way to 
the horizon by a given bulk spin network (see however [45] for a more elaborate point of view on 
this issue). The actual definition of the black hole entropy is performed by tracing out over the 
bulk states to obtain a density matrix describing a maximal entropy mixture of surface states with 
eigenvalues in the area interval of our choice. 

In order to obtain the statistical entropy in the area microcanonical ensemble the relevant com- 
binatorial problem in the U{\) framework, as explained at length in [12], consists on counting the 

sequences (cj) of non-zero elements of Z# satisfying c\ -\ \- cn = 0, and such that the condition 

cj = —2mi (modk) is satisfied for permissible spin components (mi). In this context we say that 
such a sequence of mj labels is permissible if there exists a sequence of non-vanishing spins (_//) 
such that each mj is one of the spin components of jj and 



a k 



8 < a LQG (j/,m 7 ) = 8?n4£ VJ/U7+T) < a k . (4.2) 



In practice, the counting of c-labels amounts to the determination of the dimension of the Hilbert 
subspace of ^jf or describing the black hole degrees of freedom. 

A final, and somewhat unexpected, advantage of introducing an area interval is the possibil- 
ity of simplifying the actual computation of the entropy by restating the combinatorial counting 
problem in simpler terms. Specifically, by employing an interval of the form [O,*^], Domagala and 
Lewandowski [14] proved that the black hole entropy can be obtained by the following prescription 
involving only the bulk labels mp. 

DL-black hole entropy [14]: The entropy 5^J of a quantum horizon of classical area at according 
to the Ashtekar-Baez-Corichi-Krasnov (ABCK) framework is 

S^ t (a k ) = log£l DL (a k ), 

where QP L (a k ) is 1 plus the number of all the finite, arbitrarily long, sequences (m\, . . . ,m^) of 
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non-zero half integers, such that the following equality and inequality are satisfied: 

N 

J>/ = 0, (4.3) 

7=1 

|/KI(W + 1)<^. (4.4, 

The extra term 1 above comes from the trivial sequence. 

To show that this simplification is possible one has to build a bijection between the sequences 
(c/) satisfying the conditions made explicit above, and the sequences (mi) that appear in the pre- 
ceding prescription. By doing this we show that it is equivalent to count the sequences (cj) or 
the (mi). Let us suppose first that we are given a sequence (mi) satisfying (4.3) and (4.4). If 
we take now c/ = — 2m/(mod/c) and ji = \mi\ the sequence (cf) obviously satisfies the conditions 
c\ H \-cn = and a = — 2mi(modk). Also, the sequence (mi) is trivially permissible (just con- 
sider ji = \mi\). Conversely, if we are given a sequence (c/) satisfying the conditions necessary for 
it to be counted in the computation of the entropy we can find a unique sequence (mi) satisfying 
(4.3) and (4.4). To see this we follow [14]. Suppose that we are given one of the prequantized 

values for the area corresponding to the CS level k. Now the condition ci -\ \- c^ = implies 

that 

N kL 

Z m ' = ^> ieZ' (4-5) 

7=1 Z 

In order to find a permissible sequence of spin components (mi) the following chain of inequalities 
is useful 

N 



a k = 4nyf P k > %nyl\ £ y/MJi + K 



7=1 

N 

> 87n4£vKI(KI+ 1 ) 

7=1 

N N 

> SKyfpY, \ m i\ > 87ry€p| £m/| 

7=1 7=1 

= A%yf P k\L\ , (4.6) 

where we have made use of (4.5) in the last step. As a consequence we see that L must vanish and, 
hence, the sought for sequence (mi) must satisfy the condition (4.3), known in the literature as the 
projection constraint. Another inequality that we can read off from the previous chain is 

\ m i\ < ~ , 

7=1 Z 

which implies that |m/| < k/2 for / = 1, . . . ,N. This restriction on the possible values for each mi 
is such that there is only a single choice for mi that satisfies the condition c/ = —2m/ (mod&) for 
the given c/. We conclude then that there is a unique permissible sequence (mi) associated with the 
given (a) and also that it must satisfy the projection constraint. 
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Notice that the entropy S^j(a) is defined only for area the prequantized area vaiues a k . However, 
foliowing the proposai by [i6, 17] we will extend the definition to arbitrary values of a G [0,°°) by 
just requiring that 

|^WTT)<^. 

This extension is partially justified by the fact that something similar is done in the standard treat- 
ments of the microcanonical ensemble and also because the detailed form of the area spectrum of 
loop quantum gravity is one of the main predictions of the formalism. In the following we will 
write areas in units of 4ny£p unless stated otherwise. It is important to stress at this point that the 
foregoing rephrasing of the combinatorial problem that must be solved to compute the entropy does 
not change the fact that what we really want to count are some Chern-Simons states on the horizon. 
The use of the mi labels is a (very) useful simplifying device but nothing more. 

A reasoning along the lines presented above leads to the statement of the combinatorial prob- 
lem in the SU (2) case. Now the entropy is obtained by counting states in the SU(2) Chern-Simons 
phase space. The quantum boundary condition and the introduction of an area interval lead to a 
combinatorial problem in which one has to count sequences of spin labels ( //), associated with the 
spin network edges piercing the horizon with a degeneracy factor related to the dimension of the 
invariant subspace of Inv(<g>/ [/'/]). Specifically: 

ENP-black hole entropy. The entropy 5^ p (a^) of a quantum horizon of the classical area a k = 
4ny£pk (when y < y/3) is defined as 

S^(a k )=lo g a^(a k ), 

where £2 ENP (a/t) is 1 plus the number of all the finite, arbitrarily long, sequences (j\, . . . ,_/#) of 
non-zero half integers ji satisfying 



N 

I 

7=1 



fccy^EVK// +i)<«* 



and counted with a multiplicity given by the dimension of the invariant subspace Inv((g>/[/7]). 

A nice feature of this approach is the fact that one can ultimately think of the spin network 
state representing a quantum black as one in which all the edges that pierce the horizon are bundled 
together in a single vertex labeled with an intertwiner belonging to the invariant subspace of (g)j [jj] ; 
(see [46] for more details). Moreover, it is important to point out that the most relevant part of 
the counting problem in this setting essentially coincides with the heuristic proposal by Ghosh and 
Mitra [47]. 

4.1 Combinatorics: precision countings 

In this section we describe step by step the resolution of the combinatorial problems that must 
be solved to compute the entropy according to the U( 1 ) prescription (a simple modification of the 
procedure can be used to solve the problem for the SU(2) case and other models [47, 48]). The 
steps that must be followed are [22] : 
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Step 1. Fix a given value of the area a and find the possible choices of half integers \mj\ ^ 
satisfying 

N a 

7=1 Z 

Notice that we are considering the possible choices of \mj\ as the elements of a multiset (and hence 
there is no ordering of the labels). In other words, at this stage we only find out how many times 
each spin component appears. 

Step 2. Count the different ways to reorder the previous multisets. 

Step 3. Count the different ways of introducing signs in the sequences (|m/|) in the previous step 
in such a way that £ 7 mj = 0. 

Step 4. Repeat this procedure for each area eigenvalue smaller than a and finally add the number 
of sequences obtained in each case. 

The solution to the combinatorial problem described above can be encoded with the help of a 
generating function. The coefficients of its expansion give the number of configurations that must 
be counted to compute the entropy. The derivation of the generating functions follows the steps 
given above. In fact, by following them it is not difficult to find out their final form. To see how 
this is done for the different models considered in the literature see [22, 49, 50]. 

The generating functions for the U(l) and SU(2) combinatorial problems are 

/ \ - 1 

/ oo oo \ 



i=l a=l 



r, _-l\2 / oo oo k' a + \ ,-*■«-! 

^ENP/ v v \ yz-z ) V V /^ Z 



,X; 

2 \ £i<£i v z-z 

These generating functions depend, in principle, on a variable z (introduced to take into account 
the projection constraint or the dimension of the invariant subspace in the SU(2) case) and an 
infinite number of variables xi associated with some square-free numbers p\ that appear when 
the eigenvalues of the area spectrum are written in the form Y.jOj^/p] with qj G N. The infinite 
sequence of pairs of positive integers (k' a ,y' a ) are the solutions to the Pell equations (one for each 
squarefree number p\) 

(k+l) 2 - Pi y 2 = \. 

For finite areas we only have a finite number of squarefrees pj and we only need to consider the 
finite set of variables associated with them. For example [22], if we restrict ourselves to areas 
a < 18, the only variables that we need to write explicitly in the generating functions are x\, x%, *3, 

X4, x 5 , x-i, x 9 , xio, xi 8 , x 2 2, x 2 j, x 88 , x n9 , xi 56 , x m and in the t/(l) case we can write: 

G (z,Xl,X2,X3,X4,X5,X-j,X9,XlQ,Xl%,X22,X27,X%%,Xn9,Xl56,Xl9%) 

= {\-{z 2 + zr 2 )x\ -( z l6 + z- 16 )x\ 2 -{z + zr l )x 2 -{z 6 + z - 6 )x\ 
-(z 8 +z- 8 )x4 -tf + zr A )x 2 -({i+z-^xl -(z 9 +z- 9 )xl 

~(Z U +Z.- I4 )xl -( Z 3+Z- 3 )X 10 -(Z W + Z- W )X 2 , -( Z 5 +Z -5) X22 
-( Z 12 +Z- 12 )x2 7 -( Z ll+ Z -H)x 88 -( Z 13 +Z- 13 )x 119 -(z 15 +Z- 15 )x 156 -(z 17 +Z- 17 )x 198 )- 1 
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The coefficients [z°][xj' •• ■x q [ l ^}G Dh tell us the values of the number of configurations obtained 

after following the first three steps in the previous procedure for areas of the form a = q\ \/2 -\ h 

^i98\/323 (where some of the integer coefficients q\ are allowed to be zero). For example, for 
a = Ay/2 + \/3 + \/l5 we get [z°] [xfax^G^ = 24 . 

These results can be directly used to plot the exact values of the entropy as a function of the 
area and study the detailed step structure found by the authors of [16]. An example can be seen in 
Figure 2. 
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Figure 2: Plot of the black hole statistical entropy Sj^j in the £7(1) case for a range of area values between 
a = 14 and a = 18 (in units of Aityip, where 7 is the Immirzi parameter and ip the Planck length). Notice 
the staircase structure. 

The final step in the preceding list consists in adding up the numbers obtained in steps 1-3 for 
all the possible area eigenvalues. The most practical way to do this is by using Laplace transforms. 
The reason why this is so [49] is because staircase functions such as the statistical entropy S^ t 
can be written as sums of Heaviside step functions of the form 6 (a — ciq) (with jumps at the points 
a = ao) with coefficients appropriately chosen to describe the height of the steps. As the location 
and amplitude of the steps is encoded in the generating functions given before and the Laplace 
transform of a step function is simply 



3f(e(a-a );s) 



,-a s 



it is straightforward to find an integral representation for the entropy as an inverse Laplace-Fourier 
transform (the Fourier part is related to the implementation of the projection constraint in the U(l) 
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case or the counting associated with the dimension of the invariant subspace in the SU(2) one). 
The relevant expressions for the entropy are 

1 fill fXQ + ica , °° , v_l 

expS?t!r t (fl) = t^t" / / s- l {l-2Ye- s v k{k+i) cos(Dk) e as dsd(0, (4.7) 

(2ll) l l JO J.X -ioo \ jp x J 

expS^a) = * I" f^- 1 sin 2 CO (l - £ Sin( * + *> V ^^T^dO) . (4.8) 
F stat v ' {2n) 2 iJo A-o-Zoo V £j sino / 

4.2 Asymptotics 

Two relevant questions regarding the statistical entropy are its asymptotic behavior as a func- 
tion of the area (or, in other words, check that the Bekenstein-Hawking law is indeed satisfied) 
and the study of the persistence of the equally spaced staircase structure for large areas. In any 
case we want to mention at this point that the (necessary) consideration of the thermodynamic limit 
provides an important new perspective regarding the latest issue. 

The asymptotic behavior for large areas was first derived by Meissner by looking at the singu- 
larity structure of the Laplace transform (4.7) that expresses the entropy [15]. In the SU(2) case a 
similar approach [50] can be used. The final result of these analyses are the following asymptotic 



expansions for the entropy as a function of the area (measured in units of 4n£j,y) 



Sf at {a) ~ a^ h a - - log a + o(loga) , 

SZ(") ~ ara-lloga + o(loga). 

The values of a^ L and a™ p are given by the position of the largest real pole of the integrands 
for CO = in (4.7) and (4.8). By choosing the value of y in an appropriate way it is possible to 
reproduce the Bekenstein-Hawking law with the right 1/4 coefficient. It is important to realize 
that the coefficients of the subdominant corrections are independent of the value of y and differ for 
the U{\) and St/ (2) cases. These subdominant corrections are important in order to compare our 
results with the ones obtained in other approaches (for example in string inspired models). A good 
discussion of this issue appears in the interesting paper by Carlip [51]. 

As far as the issue of the persistence of the equally spaced staircase structure is concerned, the 
best approach is to partition the set of black hole configurations in such a way that the individual 
steps are singled out. Actually there is a very efficient way to do this by using generating functions 
that can be directly obtained from the ones given above [24]. In fact, it is possible to obtain gener- 
ating functions for each of the individual steps. As these steps can be understood as unnormalized 
probability distributions it is possible to find smooth approximations for them by computing the 
relevant moments (the mean and the variance). The specific form of these moments in terms of 
the integer label that enumerates the steps is such that the equal spacing between them is easily 
derived and, most importantly, the variance that measures the width of these steps grows linerly. 
This last fact implies that for large areas the steps are wider and wider so that the steps themselves 
are smoothed out in the sum (see [24] for details). Figure 3 shows a comparison of the actual value 
of the entropy and the smoothed approximation. As can be appreciated the exact and approximate 
values of the entropy match almost perfectly. 
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Figure 3: Plot of the black hole statistical entropy Sf^ t in the U(l) case for a range of area values between 
a = 14 and a = 18 (in units of Anydp, where y is the Immirzi parameter and lp the Planck length). The red 
curve shows a smoothed approximation obtained by describing the steps with error functions as discussed in 
[23]. 



4.3 The thermodynamic limit 

The necessity to consider the thermodynamic limit can be justified in two opposite looking, 
but related, ways. If, for example, one studies the statistical entropy as a function of the energy in 
a standard thermodynamical system it is obvious that it is not a sufficiently well behaved function. 
In fact, the shape of the graph corresponds to a staircase function whose derivatives are either 
zero or are not defined at all. As important thermodynamical properties of the system are given 
by derivatives of the entropy (for example, the temperature) the statistical entropy "in raw form" 
cannot be used without some kind of smoothing. If instead of looking at S stat we consider the 
partition function, defined as Y^i e ~ (where the energy eigenvalues E\ must be counted with their 
multiplicities) we face the opposite problem in the sense that, being an analytic function of /3, it 
cannot be used to study such relevant issues in thermodynamics as phase transitions. The solution 
to both problems is the same: consider the so called thermodynamic limit in which the size of the 
system (as measured by some extensive parameter like the number or constituents or the volume) is 
taken to infinity while keeping the intensive parameters of the system fixed. Deep theorems going 
back to Griffiths [44] show that the entropy in this limit has the desired features. It is a smooth 
function almost everywhere and has the right concavity properties (for reasonable interactions) to 
guarantee the stability of the system in question. Of course there are many possible and ad hoc 
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ways to smooth out a staircase function. In many elementary texts this is simply done by taking all 
the variables that appear in concrete expressions of the entropy, including the number of particles, 
as continuous. However this is hardly a unique prescription and very few general statements can 
be derived from it. 

An important difference in the behavior of the statistical and smoothed entropy can be seen in 
the subdominant corrections for large energies. This can be seen in illustrative examples such as 
the Einstein crystal that we briefly discuss now. Let us consider N distinguishable noninteracting 
harmonic oscillators of angular frequency CO with ("normal ordered") quantum Hamiltonian 

U(N)=H 1 + ---+H N , Hi = (Oa] ai . 

In the microcanonical ensemble, the number of eigenvectors of H(A^) with energy E is just the num- 
ber of (ordered) sequences (tti,.. . ,n#) satisfying Y.u=i n i = E/co. This number can be explicitly 
written in terms of combinatorial numbers 



N 

|{(m,...,» JV ):n I -e{0}UN,£n,- = £/a>}| = < 



o HE/oogn 

«r-i + L*/-j\ 

N—l J 



Therefore, the number of microstates Sl(E,N) in an energy interval [0,E] is given by 

'N+[E/a>y 



S sm (E,N) = log£l(E,N) = log v v 

To get the thermodynamic limit, let us introduce now the energy per particle e and the entropy per 
particle Ojv(e). 

E S sm (Ne,N) 1 (N+[Ne/co 
£:=—, a N (e) := = — log 

and compute the limit 

£ CO £ 

<r(e) := lim a N (e) = -log (l + -) +log (1 + -) . 

A^°° CO £ CO 

The function cr(e) is the smoothed entropy per particle in the thermodynamic limit. If we have 
now N particles the (extensive) smooth entropy would be 

E Nco E 

S(E,N)=No(E/N) = -log{l + —)+Nlog{\ + —). 

One gets the same result with the canonical ensemble. Figure 4 shows the statistical entropy for 
several values of ./V and how they approach the entropy (per particle) in the thermodynamic limit. 

We briefly discuss now the thermodynamic limit for black holes. In the LQG framework, 
where the area ensemble plays a central role, it is possible to work either with the microcanonical 
or the canonical ensembles. In the present case it is, in fact, simpler to rely on the canonical 
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Figure 4: The plot shows that, independently of the number N of harmonic oscillators the statistical entropy 
is a staircase function whose derivatives are either zero or are not defined at all. However, the entropy in the 
thermodynamic limit is a smooth, concave, function of the energy with non-trivial derivatives 

ensemble because the partition functions can be read off directly from the integral representations 
as inverse Laplace-Fourier transforms given by (4.7), (4.8). These are 

Z Uh (a) 



1 



In 



z ENP («) = \\ 

71 JQ 



1 
1 fin 

71 Jo 



21 



k=\ 



dw 
sin 2 ft)dft) 



1 - ILi ^"V^+^sin ((*+ l)w)/sin co 



Once the partition function is known, the "free energy" -obtained by a suitable Legendre transform- 
can be immediately computed [24] and used to discuss the asymptotic behavior of the entropy as a 
function of the area. In particular we find that 

S DL (a) ~ a$ L a+-loga + 0(l) a^°o. 
S™ p (a) ~ a^ p a + 0{\), a^oo. 

In the case of the DL prescription the relevant behaviors are shown in Figure 5. As it can be readily 
seen the entropy is a concave function of the area as expected from the general theorems mentioned 
above. 

5. Final comments 

The combination of the black hole model -given by the isolated horizons- with LQG methods 
provides a definite setting to study quantum black holes. Although there are several proposals to 
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Figure 5: Plot of the black hole entropy in the thermodynamical limit in the 1/(1) case. The entropy (blue, 
solid curve in the second plot) is a smooth and concave curve, as can be seen in detail in the first plot. The 
asymptotic approximation (red, dashed curve) describes the behavior of the entropy with good accuracy. 
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carry out the quantization of these models, there is a robust formalism to study them all and extract 
the relevant physical information. In particular, the combinatorial problems that must be solved to 
compute the entropy have been thoroughly investigated and understood. The results obtained in 
all cases differ in some details (that may be relevant from a physical perspective) but some central 
predictions are robust, in particular with regard to the Bekenstein-Hawking law. 

The most important issue at this point is to take into account the dynamics. This will be 
necessary to understand, for example, Hawking radiation (for recent work on this subject see [52], 
other initial proposals appeared in [28]). We expect that the methods and ideas explained in the 
present paper will play an important role in the description of black hole evaporation. A possible 
setting where it seems possible to advance in this direction is the consideration of the models 
inspired in the Brown-Kuchaf approach to the introduction of an extrinsic time variable. These 
have been already described within the LQG framework [53, 54]. Our expectation is that the 
general ideas and methods described here will be useful in future developments in the field and may 
help extract concrete predictions in the realm of black hole physics from loop quantum gravity. 
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